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ABSTRACT 

We address the difficult question of inferring plausible node 
mobility based only on information from wireless contact 
traces. Working with mobility information allows richer pro- 
tocol simulations, particularly in dense networks, but requires 
complex set-ups to measure, whereas contact information 
is easier to measure but only allows for simplistic simula- 
tion models. In a contact trace a lot of node movement 
information is irretrievably lost so the original positions and 
velocities are in general out of reach. We propose a fast 
heuristic algorithm, inspired by dynamic force-based graph 
drawing, capable of inferring a plausible movement from any 
contact trace, and evaluate it on both synthetic and real-life 
contact traces. Our results reveal that (i) the quality of 
the inferred mobility is directly linked to the precision of 
the measured contact trace, and (ii) the simple addition of 
appropriate anticipation forces between nodes leads to an 
accurate inferred mobility. 

Categories and Subject Descriptors 

C.2.1 [Computer-Communications Networks]: Network 
Architecture and Design — Store and forward networks; 1.6. m 
[Simulation and Modeling]: Miscellaneous 

General Terms 

Algorithms, Measurement 

Keywords 

Delay Tolerant Networks, Graph Animation, Movement In- 
ference 

1. INTRODUCTION 

In the disruption-tolerant network (DTN) paradigm, mo- 
bile communication devices undergo a sequence of connec- 
tions and disconnections from other devices forming contact 
opportunities [t]. Despite the growing interest in exploiting 
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these contact opportunities for disseminating information 
under conditions when more traditional approaches are ei- 
ther impractical or unfeasible, there have been few real-life 
DTN deployments [20[ |17| . Instead, most evaluations of new 
protocols and designs have been done through simulations 
based on either synthetic mobility models or real-life contact 
traces. Neither is fully satisfactory. 

On the one hand, synthetic mobility models give full knowl- 
edge of the mobility and therefore allow for simulation of 
the specific features of radio channels (e.g., interferences and 
hidden stations) but do not accurately represent real-life 
(in particular human) mobility. On the other hand, contact 
traces are assumed to accurately represent real-life mobility 
but all geographical information is lost and simulators must 
make very simplistic assumptions on the communication 
channel (e.g., a node may only communicate with one of 
its current "contacts" at any given time |14l). A way out 
of this alternative could be the use of GPS measurements 
of human mobility [16|. Unfortunately, these are quite rare 
(only one such trace [16] on CRAWDAD [l], as opposed to 
at least 5 Bluetooth contact traces). Furthermore, they are 
often unusable as contact opportunity information because 
the distances between devices are too great. Even if large 
scale GPS measurements |2| were able to achieve a density 
allowing one to consider geographic proximity as contact 
opportunities, they would still suffer from other limitations, 
such as not working indoors. 

What if the information from the contact traces were suffi- 
cient to infer plausible node mobility? The benefits would be 
twofold. First, being able to visualize node movement is in 
itself valuable, as it confers an intuitive understanding of the 
trace dynamics that can get lost in statistics. Second, using 
the inferred movement instead of simply contacts history 
would allow for a much finer simulation of the radio channel, 
particularly for dense contact traces fisl, while retaining the 
realism captured by the contact traces. This paper makes 
the following contributions: 

• We define and discuss the problem of inferring plausible 
node mobility only from their contact information. To 
the best of our knowledge, this is the first time such a 
problem is addressed. 

• We propose a formal definition of the problem as a 
system of non-linear inequalities. 

• We describe and evaluate, both on synthetic and real- 
life contact traces, a heuristic but practical and effective 
method of inferring the movement of the nodes. 



Table 1: Some wireless contact traces. 





Duration 


Devices 


Sampling 




(days) 


(#) 


(seconds) 


MIT 


365 


100 


600 


Infocom 


3 


41 


120 1 


Rollernet 


0.125 


62 


15 


PMTR 


15 


44 


1 



When only measuring the contact opportunities from an 
experiment with mobile devices, a lot of information is irre- 
trievably lost. Consider a simple example with two nodes. 
When in contact, we can roughly locate them relatively to 
one another. However, when the time elapsed since the lat- 
est contact (inter-contact time) increases, the information 
regarding their relative distance decreases. After a while, 
it becomes difficult to say if they are still somewhat close 
or if they have gone in completely opposite directions. In 
a dense network, the higher the contact intensity, the more 
constrained our problem is. Although it is difficult to infer a 
mobility that is strongly correlated with the original mobility, 
we show in this paper that it is possible to propose a plau- 
sible mobility, i.e., one that would have generated the same 
contact trace. This is examined in more detail in Section |3l 
Since the ultimate goal is to improve simulations, inferring 
the exact mobility is not required. All we need is an inferred 
mobility that leads to better predictions. 

In the next section, we position our paper in comparison 
with prior work. In Section|3j we formally define the problem 
of inferring mobility from contact traces and discuss its 
challenges. In Sections |4] and [5] we respectively propose and 
evaluate a heuristic approach to efficiently solve our problem. 
Finally we conclude our work and discuss the path ahead in 
Section [6l 

2. RELATED WORK 

Delay/disruption-tolerant networks (DTN) 7 arise when 
lack of end-to-end connectivity, rapidly changing topology, 
and/or potentially long communication delays render tra- 
ditional mobile ad-hoc networks (MANET) approaches un- 
feasible [s]. Such networks encompass a vast spectrum of 
situations ranging from inter-planetary communications [2] 
to hop-by-hop data forwarding between portable devices to 
supplement an infrastructure for content dissemination [13| . 
In DTNs, node mobility can be exploited to increase the 
network capacity while compromising on delays by using a 
message store- and- forward paradigm instead of the usual 
packet switching. [12| . 

Opportunistic mobile networks are a class of DTNs in 
which no knowledge of the future mobility of nodes is as- 
sumed. For example, this is the case of a network formed by 
the direct contact opportunities of hand-held devices, such as 
smartphones. Contact opportunities between mobile devices 
could be used either to replace or assist a wireless infrastruc- 
ture for the dissemination of a given content. For example, 
loannidis et al. studied how to combine content pushing from 
a source in the infrastructure with opportunistic forwarding 
among subscribers in a way that ensures perceived content 
freshness from the subscribers while keeping the load on the 
infrastructure as small as possible [ 13] . 

A lot of effort has therefore gone into measuring human 



mobility, or at least the contact opportunities generated by 
human mobility. Direct measurements of human GPS coordi- 
nates have, to the best of our knowledge, only been performed 
by Rhee et al. in their work on human mobility models |16| . 
These GPS measurements give accurate and fine-grained 
(10-second period) information but unfortunately only work 
when outdoors. Measurements of human contact opportu- 
nities overcome this indoors limitation but forgo location 
information. In the Reality Mining experiment conducted at 
MIT, each participation had a special application running on 
her/his mobile phone that captured proximity information 
from 100 subjects over an academic year j^. The Haggle 
project used Intel iMotes to capture the contacts between 
participants of the Infocom 2005 conference [sj. Both of 
these experiments relied on periodic Bluetooth scans. 

While Bluetooth has the advantage of being widely avail- 
able, its scanning mechanism is too slow to effectively detect 
all contact opportunities. Indeed, the longer the sampling 
period (respectively 600 and 120 seconds for the MIT and In- 
focom traces) , the more likely temporary link failures or short 
contacts will be missed. The Rollernet experiment, which also 
uses iMotes to capture the contacts in a roUerblading tour, 
was able to bring the sampling period down to 15 seconds [18| . 
For finer measurements, a different beaconing method must 
be used. For example, Gaito et al. designed a specific device, 
a Pocket Mobility Trace Recorder (PMTR), and were able 
to measure contact opportunities every second [Tl], but the 
traces are not yet publicly available. Table [l] compares these 
different traces. As we will see in Section |4j shorter sampling 
periods translate into more constraints on our mobility in- 
ference problem, which in turn lead to solutions that better 
match the original mobility. 

Inferring node mobility based solely on contact information 
is an open problem that has not yet received any attention 
in opportunistic networks. However, some similar questions 
have been addressed in other contexts. In wireless sensors 
networks, sensors can estimate their position relatively to a 
small number of anchor nodes (typically equipped with a GPS 
receiver) using a variety of distance measurement techniques 
based on received signal strength or differences in beacon 
timings [15] . In our approach, we must also rely solely on 
contact information and assume no low-level information on 
distances between nodes. Furthermore, unlike most wireless 
sensor networks, our nodes are all mobile. Finally, and 
perhaps most importantly, being free of decentralization or 
real-time requirements, our calculations take place offline 
and with full knowledge of future contacts. 

3. INFERRING MOBILITY 
3.1 Problem definition 

Let us consider a fixed number of mobile nodes. A contact 
trace is the list of contact events that occur between these 
nodes. Each event is recorded as a quadruplet consisting 
of the identity of both nodes, the instant when the contact 
is first established, and the instant when the contact goes 
down. 

Real-life traces are noisy. In real-life traces, depending on 
the scanning period and the choice of radio technology (e.g., 
Bluetooth, ZigBee), a number of contact opportunities may 
be missed, shortened, split, or merged. For example, using 
Bluetooth, neighborhood scans typically take several seconds 



and may not detect all reachable devices due to the frequency- 
hopping nature of the protocol. By using longer sampling 
periods, it becomes difficult to detect short contacts; even 
worse, a sequence of short contacts will likely be considered as 
a single long contact. Furthermore, a Bluetooth device may 
not simultaneously scan and respond to a scan. Therefore, 
many contacts will be missed simply because of the nature 
of the underlying protocol. Other wireless technologies, such 
as the custom-made Pocket Mobility Trace Recorders 
can overcome these limitations but still have to contend 
with the traditional wireless issues such as interferences or 
hidden terminals. For these reasons, real-life traces must be 
considered noisy. 

Synthetic contact traces. Contact traces can also be 
extracted from synthetic mobility models by simulating a 
beaconing protocol or using a simple proximity-based model 
(i.e., a contact exists when two nodes are in transmission 
range of each other). Traces obtained in this fashion can be 
considered perfect, in the sense that we have full control over 
all parameters and all contact opportunities are recorded. 
We will use this approach to evaluate our heuristic algorithm 
proposed in Section [4] 

Additional topology information. Some nodes could 
have fixed and known positions, such as base stations in a 
3G of Wi-Fi network, which enables us to place other nodes 
relatively to them. This approach has been well studied 
for node positioning in wireless sensor networks where some 
sensors have GPS capabilities [l5]. When all nodes are mobile, 
relative positioning information may still be available. For 
example, in the RoUernet experiment ^18j, an iMote was given 
to a member of staff that remained at front the roUerblading 
tour, while another was given to someone who stayed at the 
back. All other nodes in the trace must therefore be placed 
between these tail and head nodes. Finally, we could suppose 
that only the initial positions of the nodes are known. For 
synthetic traces, this information is readily available. 

Plausible mobility. Since we cannot hope to recover the 
exact initial mobility from a pure contact trace, we define 
the concept of plausible mobility. In order to be plausible, 
the inferred movement must (i) realistically limit the speed 
of the nodes and (ii) possibly produce the original contact 
trace, i.e. nodes in contact must be within transmission 
range of each other while nodes not in contact should be 
beyond transmission range. 

In the end, our objective is to develop an algorithm that 
takes a contact trace and some additional information (e.g., 
fixed or relative positions) as input and generate a plausible 
movement trace as output. 

3.2 Evaluation framework 

When using synthetic mobility models, we initially have 
information of the nodes' mobility. From this, we can extract 
a contact trace, which we will use as input to our mobility 
inference algorithm. The inferred mobility can be compared 
to the original mobility, but can also be used to extract an 
inferred contact trace, which in turn can be compared to the 
original contact trace. When using real-life contact traces, we 
can no longer compare original and inferred mobility, but it 
still possible to compare the contact traces. This is summed 
up in Fig [l] 

We consider two ways of evaluating a mobility inference 



Original movement 



extract contacts 



Original contact trace 




Inferred contact trace 



Figure 1: Evaluation framework. 

method: one comparing original and inferred mobility, and 
the other comparing original and inferred contact traces. 

3.3 Formalization 

In this section, we describe what would constitute an ideal 
inference of mobility. The constraints defined below will guide 
us in the choice of the parameters for the heuristic approach 
proposed in Section [4] Since the input is a contact trace, 
complete knowledge of past, present, and future contacts is 
assumed (offline inference). 

3.3.1 Definitions 

Let A'^ be the number of mobile nodes in the contact trace. 
These nodes move on a 2D plane, have a maximum speed 
Wmax and a transmission range r. Let {xi{t),yi{t)) be the 
coordinates of node i € {1, • • • For the pair of nodes 

(i, j), let dij{t) denote the distance between i and j at time t. 
Furthermore, Tj.{t) and T^j{t) denote, respectively, the time 
at which the next contact between i and j will appear and 
the time when the current contact will end. 

At any time t and any time interval Ai, the maximum 
node speed fmax imposes the following constraint on the 
positions of any node i: 

J {xi{t + At) - x^(t)f + {yi{t + At) - yi{t)f < ''-'max 

At. 

(1) 

This constraint imposes that, given a valid solution at time 
t, and a short time interval At, the next valid position at 
time t + At should be very similar. 

3.3.2 Case 1: Synthetic contact traces 

In a synthetic contact trace, a contact appears when the 
distance between two nodes is less than r and breaks when 
it is greater than r. Since we know when current contacts 
are going to break and when new ones will appear, we can 
further constrain nodes' positions. Indeed, nodes must get 
closer to each other before the contact appears and move 
away from each other before it goes down. This is illustrated 
in Fig. |2j where the contact at time t between nodes i and 
j will go down at T^j{t) and a contact between i and k will 
appear at time Tji,{t). As t approaches T^^(t), node j must 
get close to moving out of i's transmission range. Relatively 
to i, j must be able, given its maximal speed Umax to move 
out of transmission range at exactly T^-{t). Therefore, if i 
and j are in contact, the following constraint holds: 




i's transmission range 

Lower bound for j's distance to i 

Upper bound for fc's distance to i 



Figure 2: Constraints on tlie positions of nodes j and k 
relative to node i. Nodes i and j are in contact at time t 
until Tf-{t). A contact between nodes i and k will appear at 

time T^kit)- Parameters r and Wmax denote the transmission 



range and maximum speed, respectively. See Section 3.3.2 



< d,j{t) < r. 



(2) 



A similar analysis holds for the contact between i and 
k. As t approaches node k must come closer to i's 

transmission range. Relatively to i, k must be able, given 
iimax, to come within transmission range of i. Therefore, 
while i and k are coming into contact, the following constraint 
holds: 



r < d„(i) <r + u„iax (t - r7fe(t)) 



(3) 



3.3.3 Case 2: Real contact traces 

While we know that a movement satisfying constraints , 
([2|, and ([3| exists for the synthetic contact trace (i.e., the 
original synthetic movement), it is less clear for a real-life 
trace. Indeed, as we previously discussed, a real-life contact 
trace may be quite noisy and, in particular, miss many 
contacts. While this may seem like a simple relaxation of 
our constraints, it could in fact make the system unsolvable. 
Indeed, when considering real-life traces, the inclusive (i.e., 
m-contact, Eq. [2| and the exclusive (i.e., not-m- contact, 
Eq. [3| constraints no longer have the same importance. The 
inclusive constraint, based on the presence of a contact, 
can be trusted. However, the exclusive constraint, based 
on the absence of a contact, no longer strictly means that 
the distance between two nodes must be greater than the 
transmission range r. Indeed, one could imagine a node 
quickly passing by the other nodes, moving in and out their 
transmission ranges without triggering any contact detection. 
If we strictly enforce the exclusive constraint, such movements 
may no longer be possible. 



4. HEURISTIC SOLUTION 

In this section, we propose and evaluate a simple and effi- 
cient heuristic for inferring node mobility from their contact 
traces. Note that in order for it to have broad applicability, 
it should make as few assumptions as possible on the original 
mobility. 

4.1 Dynamic graph drawing 

Our heuristic approach is inspired by works in dynamic 
graph animation, even though its objective is quite different. 
Graph animation aims at (i) producing a sequence of readable 
and aesthetically pleasing representations of graphs and (ii) 
animating the transitions between successive graph layouts in 
a way that preserves the viewers' mental map of the graph [9]. 
Sample applications include visualization of communication 
networks, social networks, and software library dependencies. 
Both goals are relevant to us. Not only do we wish to infer 
a sequence of positions for each node in the contact trace 
(i.e., a sequence of connectivity graph layouts), but real-life 
mobility intrinsically produces a sequence of connectivity 
graphs in which the transitions are easy to follow. However, 
while the function of dynamic graph animation is mostly 
aesthetic, our heuristic aims at meeting the constraints set 
out in Section [3.31 

In the context of graph animation, online means that the 
graph layout algorithm is continuously running while new 
nodes or links appear and disappear on the fly, whereas 
offline means that each successive graph is laid out sepa- 
rately. The offline method makes it difficult to preserve the 
viewer's mental map during transitions, particularly long 
ones, between successive graph layouts. The online method 
procures an illusion of continuous mobility and allows for 
easy control of nodes' speed but is, in itself, insufflcient in 
our case. Indeed, when a contact occurs between two nodes, 
they may, at that time, be very far away from each other 
in the online graph animation. This leads to a link in the 
connectivity graph that will, at least temporarily, straddle 
several connected components, which cannot constitute a 
satisfactory inferred movement. 

Of particular interest to us are the force-based layout 
algorithms [lO], in which attractive and repulsive forces are 
applied to nodes according to the connectivity graph. As 
in a real physical system, the nodes then converge to a 
minimum stress (or energy) position. Force-based algorithms 
are particularly well suited to our problem because each pair 
of nodes that are in contact will tend to be geographically 
close to each other. 

Our heuristic for inferring a plausible mobility from a 
contact trace will consist of running an online force-based 
dynamic graph layout algorithm, built from the forces and 
refinements described in the next two sections. 

4.2 Forces 

As in a physical system, each node has a position, a speed, 
and an acceleration. All nodes have the same mass. Each 
node is subject to the four forces described below. The first 
three are classic, while the last is novel and necessary for the 
quality of our inference. We first present all of them before 
discussing how to set their parameters in the next section. 
Fig. [3] shows how these forces add up. Hereafter, Uij is the 
unitary vector directed from i to j. 

Attraction. Let k be one of node i's contact (i.e., a neighbor 



k 



Contacts: 

Current 

Future 

None 



Figure 3: Forces applied to node i from nearby nodes j, k, I. 
Fjj is the attractive force of j on i, FJ^ is a repulsive force, 
and F^"* is an anticipation force. F is the resulting force 
applied to i, bringing i closer to k before the {i, k} contact 




appears. Full explanation in Section 4.2 



in the connectivity graph), j attracts i with a spring-like 
force: 



Ffc, = K{d,k - lo)uik 



(4) 



where dn^ is the distance between i and k, lo is the spring's 
equilibrium position, and _fs' is a rigidity constant. This 
force contributes to keeping a pair of nodes in contact within 
transmission range of each other (right part of constraint ([2|). 

Repulsion. Each node j i pushes node i back with a 
coulomb-like force: 



if dij *C dmax-j 
if dij ^ dmax-j 



(5) 



where G is an intensity constant, eo is small strictly positive 
constant to keep this force bounded, dmax a cutoff distance 
beyond which this force no longer acts, and a a parameter 
that determines how this force's intensity decreases with 
distance. This force contributes to keeping nodes that are not 
in contact away from each other (left part of constraint (|3|). 

Drag. In order to prevent excessive speeds (cf. constraint |T])) 
and isolated nodes from moving too far away, each node i is 
subject to a drag force: 



-£>Vi, 



(6) 



where D determines how strong the drag is and is the 
current speed of node i. 

Anticipated attraction. Since we have access to the entire 
contact trace, we know which nodes will meet node i in 
the future. The idea here is to progressively make node i 
move towards its future contacts, so that when the contact 
does appear in the trace, both nodes will be roughly in 
transmission range of each other (right part of constraint ([3|). 
If node k meets i at time tik, then, at time t < tik, k attracts 
i with a spring-like force: 



K{dik - lo)e 



Uife, 



(7) 



which becomes equal, at t = tu, to the regular repulsive 
force depicted in Eq. |5] The r parameter characterizes when 
future contacts begin to have a noticeable influence on i's 
movement. 



4.3 Issues and usage 

When animating the graph, one must keep in mind that 
the goal is for each node to be within transmission range r 
of its current contacts, and outside of range of all the other 
nodes. This should be encouraged but not strictly enforced, 
as it may otherwise lead to impossible configurations for real- 
life contact traces. Allowing the possibility of inexactitudes 
(i.e., adding or removing links in the inferred contact trace) 
is the price to pay for being able to infer movement from 
real-life traces. 

A cluster of nodes can collectively have a strong repulsive 
force. As such, if another node comes into contact with one 
node in the cluster, they may never come into transmission 
range of each other, despite the attractive force. Setting a 
strong rigidity constant and setting the equilibrium length 
of the spring force to a point within the transmission range 
shown in Eq. |4] offsets this. In the rest of this paper, we use 
K — 30 and Zo = §■ In Eq. Isl we set Q = | and eo = 1. This 
ensures that the repulsive force (i) is bounded by G and (ii) 
does not decrease too quickly. Finally, we choose G in Eq. [S] 
so that the distance between two nodes in contact, absent 
all other forces, converges towards |r. 

An issue not usually addressed in the graph drawing com- 
munity is how to deal with disconnected components. Since 
we are handling DTN contact traces, we cannot avoid this 
problem, as the connectivity graph is almost always split 
into several disconnected components and many isolated 
nodes. Freivalds et al. propose laying out each connected 
component separately and then using a packing algorithm 
to place them relatively to one another 8 . However, this 
completely ignores that, in our case, the relative placement of 
connected components should not be arbitrary. Fortunately, 
our anticipated attraction force circumvents this problem by 
creating attractive forces between disconnected components 
and thus guiding their relative movement, orientation, splits, 
and merges. The value of r in Eq. [7] is a tradeoff. Small 
values of r mean that only the very short term future is 
considered for animating the contact trace, while large values 
can create so many constraints that no movement is possible. 
Good values are linked to the characteristic evolution time 
of the connectivity graph. Finally, cutting off the attractive 
force eliminates long range interactions between nodes that 
could interfere with the initially weak anticipated attraction 
forces. In the rest of this paper, we use dmax = 3r. 

5. EVALUATION 

In this section, the heuristic described above is applied both 
on a synthetic mobility model (Random Waypoint - RWP) 
where the contact trace is considered perfect and on the real- 
life contact trace with the shortest available sampling period, 
Rollernet Is]. Note that our heuristic makes absolutely no 
assumptions about the underlying mobility model. In fact, 
it is particularly poorly adapted to Random Waypoint, since 
in it the nodes try to avoid each other, whereas in RWP 
nodes pay no attention to each other. Nevertheless, we still 
manage to infer plausible movements. 

5.1 Synthetic movement 

The synthetic mobility scenario considered in this section 
consists of 50 nodes moving according to the RWP model 
on a 1000m x 1000m torus, so as to avoid any border effects. 
RWP has some well known shortcomings, such as a usually 
non-uniform steady state spatial distribution of nodes and 
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(a) Original mobility (RWP) (b) Inferred mobility 

iMl Anchor node Missing contact in the inferred mobility 
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Figure 4: Mobility snapshots of the original Random Waypoint movement (a) and the movement inferred from its contact 
trace (b) after 200 seconds. Nodes 0, 1, 2, and 3 are fixed anchor nodes. The other 46 nodes move on a torus. 



gradual speed decay |T^. However these are irrelevant in this 
work since the steady state spatial distribution on a torus is 
in fact uniform, and we only run short RWP simulations (300 
seconds) aimed solely at producing random contact traces 
with which to study our mobility inference algorithm. Fur- 
thermore, we consider that 4 nodes, hereafter called anchor 
nodes, are immobile and symmetrically distributed on the 
torus at positions (250m,250m), (250m,750m), (750m,250m), 
and (750m, 750m). These anchor nodes serve as reference 
points for placing all the other nodes. In a real-world sce- 
nario they could wireless access points or 3G base stations. 
Unless otherwise noted, nodes move at speeds chosen uni- 
formly between Im/s and lOm/s with no pause time, and 
their transmission range is 100m. For each run of this mobil- 
ity scenario, a contact trace was extracted from successive 
snapshots of the nodes' positions with a certain sampling 
period, by simply considering that any pair of nodes within 
transmission range of one another are in contact. A snapshot 
of this mobility is seen in Fig. |4^. 

Through some experimentation, the r parameter of the 
anticipated attraction force was set to 5. This means that 
future contacts start significantly pushing their nodes closer 
to each other 5 seconds before the contact actually appears. 
Smaller values meant that contacts scheduled to appear in 
the original trace would only show up later in the inferred 
trace, because of the delay until both nodes would get into 
transmission range. Greater values meant that each node 
would be attracted to a larger subset of the other nodes. 
If too big, this can result in preventing most movement. 
Additionally, the initial positions of all nodes were known to 
our heuristic. 

We first compare the inferred mobility to the original. 
While we do not expect to be able to infer the exact node po- 
sitions, the relative distances between pairs of nodes should 
be correlated. Indeed, Fig. [5] examines the correlation be- 
tween the pairwise distances in the original and inferred 
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Figure 5: Correlation between pairwise distances in a syn- 
thetic mobility trace and its inferred mobility. Correlation 
values are given in parentheses. 



mobility. Every time step, the distance between each pair 
of nodes was measured in both the original and inferred 
mobility. Fig. [5a| represents the correlation scatter plot in 
the absence of anchor nodes, while in Fig. [5b] the four anchor 
nodes are present. The presence of a few fbced anchor nodes 
clearly helps positioning the other nodes. 

One artifact of our heuristic is visible for small distances. 
Indeed, all distances below the transmission range (100m) 
in the original mobility are roughly mapped to 75m in the 
inferred mobility. This is a result of the equilibrium distance 
of two nodes in contact as discussed in Section 14.31 An 
example is visible on Fig. [4] where nodes 1, 35, and 36 are 
closely grouped in the original mobility (Fig. Hk) but appear 
as a triangle in the inferred mobility (Fig. |4p) . 

We then compare the inferred contact trace to the original. 
Every time step, we check the constraints Q and ([3|. If 
two nodes that should be in contact are not within transmis- 
sion range of each other, that contact is considered missed. 
Conversely if two nodes not in contact are within transmis- 
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Figure 6: Influence of contact sampling period on quality of 
mobility inference. Values expressed as percentage of existing 
contacts. 
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Figure 8: Inter-contact time distribution for original and 
inferred mobility in RoUernet 
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Figure 7: Influence of maximum node speed on quality of 
mobility inference. Values expressed as percentage of existing 
contacts. 



sion range of each other, an added contact is counted. In 
Figs. |6] and [7j missed and added contacts are expressed as a 
percentage of the number of existing contacts at that time. 

This is a rather strict way of comparing two contact traces, 
as even slight delays in the start or end of a contact will 
register as respectively a missed or added contact for that 
time period. Fig. [6] shows the effect of the sampling period of 
the original contact trace on the quality of our inference. The 
proportion of both added and missed contacts increases with 
the sampling period. This is due to several reasons. Firstly, 
as the sampling period increases, it becomes more difficult 
to assume that a contact present in one period but not the 
next lasted the whole period. Lacking other information, 
our heuristic does however make this assumption. During a 
time period, a given node, when pushed by the anticipated 
attraction force towards its next contact, is still restrained 
by the attractive forces of the nodes that were in contact 
with it at the beginning of the period. For longer sampling 
periods, we may overestimate the duration of many contacts 
and therefore prevent a node from moving towards its future 
contacts. This translates into both missed contacts, from 
not getting within transmission range of new contacts in 
time, and added contacts from remaining within transmission 
range of old contacts. Secondly, smaller sampling periods 
also catches short contacts that would otherwise be ignored. 
These provide many extra contacts that a node must meet 
on the way to meeting its next contacts according to the 
longer sampling period, and thus enable a much smoother 
and progressive mobility inference. 



High-frequency contact measurements are therefore neces- 
sary in order to achieve good mobility inference. This likely 
holds for any real or synthetic mobility, and any inference 
algorithm. Unfortunately, no such high-frequency contact 
traces are publicly available. 

Fig.[7|examines how the dynamic of the underlying mobility 
impacts the quality of the inference. The sampling period 
here is 1 second, and we use different values for the maximum 
speed in the Random Waypoint model. Both added and 
missed contacts increase with the maximum node speed. 
The more dynamic the movement, the harder the inference. 
Again, we expect this to hold for any combination of mobility 
models and inference algorithms. 

5.2 Real-life traces 

As discussed previously, the longer the sampling period in 
a contact trace, the more difficult it becomes to claim that 
a contact translates into a link that lasts roughly as long 
as the sampling period. Several contact traces (see Table [l| 
were considered, but in this section we focus on the RoUernet 
trace [iS], collected during a roUerblading tour around Paris, 
which due to its low sampling period (15 seconds), is the 
one that comes closest to capturing the evolution of the 
connectivity graph. 

Compared to the synthetic mobility in the previous sec- 
tion, this contact trace results from a real-life experiment and 
thus must be considered lossy. Furthermore, the nodes are 
highly mobile and the average contact duration (26 seconds) 
is barely longer than the sampling period (15 seconds). A 
look at Fig. |6] may lead one to assume that a 15-second 
sampling period is too long for a satisfactory mobility infer- 
ence. However, unlike in the previous section, these sampling 
periods are not synchronized across devices. The contact 
trace thus provides us with a continuous, though incomplete, 
stream of contact events to guide our heuristic. 

When inferring mobility for the RoUernet contact trace, 
we used the same parameters as in the previous section. 
There are no flxed anchor nodes in the trace, but the head 
and tail nodes of the roUerblading tour are known. No 
participant skated ahead (resp. behind) the head (resp. tail) 
node. The head and tail nodes were constrained to moving 
along a horizontal axis. If ever a node wants to pass the 
head node for example, then the head node is moved to 
ensure that it remains ahead. This enables all nodes to be 
placed relatively to the head and tail nodes. Due to the 
rapid contact process in the RoUernet trace, the anticipated 
attraction forces are sufScient to keep the roUerblading tour 



compact and naturally lead to the emergence of the accordion 
phenomenon 
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because the head and tail node get closer 
when the contact density increases and pushed apart when 
it decreases. The resulting mobility is aesthetically pleasing 
and helps guide intuition when working on the dataset. 

Fig.jsjplots the inter-contact time distribution for the first 
5000 seconds of both the original and inferred RoUernet con- 
tact traces. As usual, both distributions follow a truncated 
power law [s]. However, the inferred contact trace contains 
many more contacts than the experimental one. Indeed, on 
average, it adds about 36% extra contacts and misses about 
9% of the measured ones. This may seem like a lot, but 
unlike in the previous section, some of those added contacts 
may in fact be real contact opportunities that the Bluetooth 
devices simply failed to pick up. Are we introducing bo- 
gus contact opportunities? Or are we identifying contacts 
that the experiment missed? Or both? The lack of reliable 
high-frequency contact traces make it impossible to decide. 



6. CONCLUSION AND FURTHER WORK 

In this paper, we examined a new and interesting problem, 
the inference of a plausible mobility from a wireless contact 
trace. Indeed, mobility is more difficult to measure but 
enables better simulations, particularly in dense networks, 
whereas contact information is easier to measure but only al- 
lows for simplistic simulation models. Our heuristic solution, 
based on ideas from dynamic graph drawing, can animate 
any wireless contact trace, while making practically no as- 
sumptions on the mobility that produced the contact trace. 
Our results highlight the need for reliable high-frequency 
contact traces in order to extract a plausible mobility from a 
wireless contact trace. 

This work, a first approach of a difficult problem, can be 
pursued in several directions. A more comprehensive study 
of our proposed heuristic must be undertaken using different 
mobility models and scenarios. In particular, the time to 
reach a threshold of accuracy when starting without knowl- 
edge of initial positions in the experiments of Section. |5.1| 
should be examined. Furthermore, other heuristics, perhaps 
closer to the constraints in Section [3. 3| should be explored. 
Finally, the idea that infering mobility from a contact trace 
before running simulations on it ought to lead to more re- 
alistic results than just using the contact trace alone, must 
be properly tested. For example, one could compare the 
performance of various network protocols in an opportunistic 
network based on simulations, with complete knowledge of 
the mobility, with only contact information, and with the 
plausible mobility infered from the the contacts. The results 
on the infered mobility should be closer to those using the real 
mobility than to those using only the contact information. 
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